Inverse kinetic theory for quantum hydrodynamic equations 
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We propose a solution for the inverse kinetic theory for quantum hydrodynamic equations as- 
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I. INTRODUCTION 

A striking feature of standard quantum mechanics 
(SQM) is its analogy with classical fluid dynamics. In 
fact, it is well-known that the Schrodinger equation 0| 
is equivalent to a closed set of partial differential equa- 
tions for suitable real-valued functions of position and 
time (denoted as quantum fluid fields) . The quantum 
hydrodynamic description obtained in this way, which 
has been studied by several authors 0, IS S IS S S 
0,0,0,01, affords a straightforward physical interpre- 
tation in terms of a classical fluid. This is obtained by 
identifying the classical fluid fields (number density and 
fluid velocity) respectively with the quantum probability 
density in configuration space and the quantum proba- 
bility current (or quantum fluid velocity). In particular, 
the quantum hydrodynamic equations can be viewed as 
the equations of a classical compressible and non- viscous 
fluid, endowed with potential velocity and quantized ve- 
locity circulation. 

The analogy (with classical fluid dynamics) has moti- 
vated in the past efforts to formulate phase-space tech- 
niques based on various statistical models of quantum 
hydrodynamic equations. These works, although based 
on different approaches, share the common view that the 
quantum state corresponds to an underlying ("subquan- 
tum") statistical description of some sort. 

Starting from the pioneering work of Wigner 0, 0| , 
phase-space techniques typically require that the quan- 
tum fluid fields or the wave-function itself be repre- 
sented in terms of, or associated with, appropriate phase- 
space functions. These are usually identified with quasi- 
probabilities although formulations based on pure 
probability densities are also available (see Ref . [l^ . sec. 
IV. 3). The procedure of representing quantum states 
by quasi-probabilities in phase space is closely related 
to the phase space formulation of SQM mechanics based 
on the noncommutative product known as Moyal product 
[lol ll^. However, quasi-probabilities often become sin- 
gular or negative in the so-called "full quantum regime" 
|l8l |. i.e., when the potential becomes nonlinear. This fact 
has motivated in the past the search of alternative phase- 



space representations of the quantum state. Among such 
approaches, we recall the phase space representation of 
SQM due to Torres- Vega and Frederick [3, , in which 
the wavefunction is extended to phase space, V'(r), 
is assumed to obey an appropriate " Schrodinger equation 
in phase space" . 

Another class of statistical approaches is represented 
by the attempt to interpret SQM in the framework of 
an underlying subquantum statistical process. An ex- 
ample is p rovided by stochastic models (see for example 
Refs. HH m m Hi 113 and the review paper H^), in 
which the underlying particle dynamics is governed by 
stochastic differential equations, such as those describ- 
ing the nondissipative quantum Brownian motion, which 
traditionally is described by Fokker-Planck or diffusion 
equations. Such equations, generally, lead again to quasi- 
probability distributions which permit to "reconstruct" 
the Schrodinger equation only approximately (namely in 
an asymptotic sense) and under "ad hoc" initial condi- 
tions, since quasi-probability functions may become, in 
general, invalid for non-Markovian processes with arbi- 
trary noise correlation. Although extensions of the clas- 
sical theory of Brownian to quantum domain have been 
proposed which permit the formulation in terms of a 
well-behaved true probabilistic description , the prob- 
lem of these approaches remains that of explaining the 
origin of such stochasticity H^, which is precisely 
their weakness. 

A second type of statistical approaches is provided by 
classical subquantum kinetic models in which the un- 
derlying particle dynamics is purely classical, i.e., the 
Schrodinger equation is assumed to be based on classi- 
cal kinetic theory. An example of this type is provided 
by the approach due to Kaniadakis 30, 31], in which 
each quantum particle is assumed to be composed by 
N ^ 1 identical subquantum interacting classical par- 
ticles (monads), to be described by means of a classical 
kinetic equation. Also in this case the Schrodinger equa- 
tion is - at best - recovered only in an approximate sense. 
Several aspects of this formulation, however, remain un- 
clear, which include - among others - the problem of the 
closure of the moment equations, the specification of the 



2 



initial and boundary conditions for the kinetic distribu- 
tion function and the conditions of convergence to the 
Schrodinger equation. 

In this work we intend to investigate the problem of 
searching an inverse kinetic theory for the Schrodinger 
equation, i.e., a kinetic theory yielding exactly, by means 
of a closed set of moment equations, the quantum hydro- 
dynamic equations. By definition, such a theory should 
be non-asymptotic and able to yield the correct equations 
independently of the initial and boundary conditions for 
the wavefunction, to be assumed suitably regular. Fur- 
thermore, it should hold also for arbitrary, but suitably 
smooth, initial conditions for the kinetic probability den- 
sity and the moment equations should form a complete 
system of equations, namely satisfy a closure condition. 

The construction of an inverse kinetic theory for the 
Schrodinger equation involves possibly also the identifi- 
cation of an underlying classical dynamical system, in 
terms of which all relevant observables and correspond- 
ing expectation values can in principle be advanced in 
time. 

From the mathematical viewpoint the problem can be 
set, in principle, for arbitrary fluid equations, an exam- 
ple being provided by the inverse kinetic theory recently 
developed for the incompressible Navier-Stokes equation 
by EUero and Tessarotto HE 1^ (hereon referred to 
as ET). A basic prerequisite for the formulation of an in- 
verse kinetic theory of this type is, however, the proper 
definition of the relevant quantum fluid fields and their 
identification with suitable momenta (to be denoted as 
kinetic fluid fields)^ which include the kinetic tempera- 
ture as well as the related definition of directional tem- 
peratures (see below). 

In this sense, the case of Schrodinger equation is pecu- 
liar because, as is well-known, its related fluid equations 
apparently depend only on two quantum fluid fields, re- 
spectively, to be identified with the observables quantum 
probability density and the quantum fluid velocity, while 
the notions of quantum temperature and directional tem- 
peratures (to be identified with the corresponding kinetic 
moments) remain in principle arbitrary. In this paper we 
intend to propose a possible definition of these observ- 
ables which is related to the kinetic interpretation of the 
Heisenberg theorem. 

The problem is not merely of interest for theoreti- 
cal and mathematical research, but has potential rel- 
evance also for the understanding of the fiuid descrip- 
tion of quantum mechanics and of the underlying statis- 
tical models. Our motivation, in particular, is to ex- 
ploit the analogy between classical and quantum sta- 
tistical mechanics, in order to prove that the quantum 
observables and the fluid flelds can formally be repre- 
sented by means of a purely classical statistical model. 
Although the mathematical equivalence should not too 
hastily be regarded as implying physical equivalence of 
the two formulations, this suggests that some relevant 



classical reasonings can be transferred to SQM. This in- 
volves, in particular, the formal description of SQM by 
means of a classical dynamical system (to be denoted 
as phase-space Schrodinger dynamical system) which de- 
scribes the dynamics of a set of virtual subquantum par- 
ticles interacting with each other only by means of an 
appropriate mean-field interaction. 

Such a classical description is realized by means of 
an appropriate form of the correspondence principle 
(denoted kinetic), whereby the physical observables, 
quantum fluid fields and quantum hydrodynamic equa- 
tions are respectively identified with appropriate ordi- 
nary phase-space functions, velocity moments of the 
(subquantum) kinetic distribution function and appro- 
priate moment equations obtained from the inverse ki- 
netic equation. 

Here we intend to show, however, that in principle in- 
finite solutions to this problem exist, namely the inverse 
kinetic theory is intrinsically non-unique. Such a feature 
is not surprising since these kinetic models may be viewed 
simply as examples of the infinite admissible, and phys- 
ically equivalent, mathematical descriptions of physical 
reality. However, despite such a non-uniqueness feature, 
it turns out that the functional form of the kinetic equa- 
tion and the mean-field force that defines the streaming 
operator, by suitable prescriptions, can both be uniquely 
determined. 

Another implication of the theory is to show in a sim- 
ple way that the Bohmian program, of reproducing the 
predictions of SQM within a framework in which parti- 
cles have definite Lagrangian trajectories ^IHI^t has by 
no means a unique solution, since there are infinite equiv- 
alent possible realizations of the Schrodinger dynamical 
system which satisfy the quantum hydrodynamic equa- 
tions. As is well-known 11, 12J, this means that there 
are in principle infinite alternative theories to Bohmian 
Mechanics, which attribute deterministic trajectories to 
the particle wave-function and are equivalent to SQM. 

The paper is organized as follows. 

The mathematical setting of the hydrodynamic de- 
scription of SQM is recalled in Sees. II-IV. In particular, 
the quantum hydrodynamic equations for the quantum 
fiuid fields {/, V} are posed in Sec. II and III, while 
the definition of the remaining quantum fiuid fields, pro- 
vided by the quantum directional temperatures (Tgjv/.i, 
for i = 1, 2, 3), is given in Sec. IV, together with their rela- 
tionship with the Heisenberg inequalities. The construc- 
tion of the inverse kinetic theory is treated in Secs.V and 
VI. In particular, the general form of the inverse kinetic 
equation is laid in Sec.V, while in Sec. VI the mean-field 
force is determined explicitly. The main results of theory 
are summarized in three theorems which determine the 
explicit form of the inverse kinetic equation and its basic 
properties. 
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II. HYDRODYNAMIC DESCRIPTION OF NRQM 



In this section we intend to recall the well-known 
fluid description of non-relativistic quantum mechanics 
(NRQM) , based on the property of the Schrodinger equa- 
tion to be equivalent to a complete set of fluid equations. 
For the sake of clarity let us introduce the basic defini- 
tions and the mathematical formulation of the problem. 

In this paper we shall consider, in particular, the case 
of a system of spinless scalar particles (bosons) described 
by a single scalar wavefunction '0(r,t),with associated 
probability density 



f = mr,t)\' 



(1) 



requiring that both are defined and continuous in 51 x / 
and in addition that / is strictly positive in fi, il denot- 
ing the closure of il. Hence, by assumption, / can only 
vanish on the boundary SQ (i.e., in the nodes r„ G 6fl 
where /(r„,t) = 0) and must satisfy the normalization 



dr/(r,t) = l. 



(2) 



In addition, / and ip are respectively assumed single- 
valued and possibly multi- valued in x /, with ?/> at 
least of class C^^^^-'(ri x /). For definiteness, we shall also 
assume, without loss of generality, that SI is a connected 
subset of M''^ and ^p{r, t) belongs to the functional space 
{■0} , to be identified with the Hilbert space of complex- 
valued functions which arc square-integrable in f2. The 
TV— body wave-function i/'(r, i) is required to satisfy in 
the open set ft x I the vector Schrodinger equation 



d 



(3) 



where H — Hq + U is the iV— body Hamiltonian operator. 
Here, denoting V(,) = = - J2 ^^^7 and U 

j=l.N 

are respectively the free-particle Hamiltonian (kinetic en- 
ergy) and the interaction potential, to be identified with 
a real function defined and suitably smooth inQx I. For 
well-posedness, appropriate initial and boundary condi- 
tions must be imposed on ip{r,t). The initial conditions 
are obtained by imposing for all r e 



-0(r,io) = -0o(r), 



(4) 



where ipo is a suitably smooth complex-valued function. 
To specify the boundary conditions, we first notice the 
boundary set Sn can always be considered prescribed. 
The boundary conditions can be specified by imposing 
Dirichlet boundary conditions on S^l. This requires Vr^ G 

sn 



(5) 



lim V(r,t) = V,„(r,,i), 

r— >ra 

where V(r, t) is the quantum velocity field 

h 



(6) 



V(r,i) 



2mi |'0(r, t)\ 



{V'*V?/' - i/'VV'*} . (7) 



Here the complex function il>w{v^^t) and the real vector 
function V^, (r^ , t) are prescribed, suitably smooth func- 
tions. To specify the value of /(r, t) on (5S1, let us require 
that there results additionally 



/ drV/(r,t) = 0. 

JQ. 



(8) 



In all such cases Eq. implies that there must be Vr^ e 
50. 



/(r,,i) = |V™(r,,t)|' = /„>0, 



(9) 



where /o is either a constant, whose value may still de- 
pend on the specific subset, or at most is a function 
fo{t) to be assumed suitably smooth Mt G /. Hence, 
the points of 50. are not necessarily nodes. However, 
if is an improper point of M?"^ (hence, fl is assumed 
to be an unbounded subset of R'^^), since it must be 



lim I 



J(r,t) 



0, is necessarily a node, i.e., 

fo = 0. (10) 



This implies for consistency also 



|r| 



hm i/;^(r, t) — 0. 



(11) 



The set of equations JSJ , 0) , © , ® together with © or 
(|10|l and Hll|l , defines the initial-boundary value problem 
for the Schrodinger equation [SE problem). The solution 
of the SE problem, t/;, must be determined in an appro- 
priate functional space, to be suitably defined (see for 
example R,ef.|8^). 

The set of hydrodynamic equations corresponding to 
the Schrodinger equation are well-known 2, _3j .4, 8] and 
follow immediately from the exponential representation 
(known as Madelung transformation 0) 



e R . 



(12) 



where {/, S*} , denoted as quantum fluid fields, are re- 
spectively the quantum probability density and the quan- 
tum phase- function (also denoted as Hamilton-Madelung 
principal function). Eq. 1)1211 is manifestly defined only 
in the set in which results / > (i.e., in the configuration 
space O). We stress that in principle S{r,t) remains "a 
priori" unspecified on the subset the boundary SO. where 
/ = (subset of the nodes r„). This indeterminacy, how- 
ever, is eliminated by requiring that everywhere in 50, 
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S'(r, t) can be prolonged on the same set by imposing 
Vr„ G Sfl : 



S{rn,t) = lim S{r,t). 



(13) 



Hence, the real functions {/, 5*} can both be assumed 
continuous in x / and at least C'^'^'^-' (fix I). Obviously, 
S{r,t) is defined up to an additive constant 27rfc, being 
fc G Z, while S itself is generally not single- valued. In 
addition, if ■0 is single- valued, it is obvious that S must 
satisfy a well-defined condition of multi-valuedness. In 
fact, in this case on any regular closed curve C of il, S 
it must result 



(14) 



The initial conditions to be satisfied by the quantum 
fiuid fields {/, S} , stem from Eq. Q and read: 



/(r,io) = /„(r). 



S{r,to) = So{r) morf(27r). 



(20) 



(21) 



Instead, the boundary conditions implied by Eqs. 
(El),®,©,© together with © or and read 
respectively Vr^ G Sil 



f{'rs,t) = fw{rs,t) 



(22) 



where n is an appropriate relative number |35| . Intro- 
ducing the single-valued potential velocity field, defined 
in fl X I, 



1 



V(r,i) = — VS'(r,t) 



(15) 



this yields the well-known condition of quantization of 
the velocity circulation 



K = J dl-V{r,t) = 



2'Knh 



m 



(16) 



Hence, by denoting -gj = ^ + V • V (convective deriva- 
tive) , it follows that in the open domain Ox/ (where by 
definition / > 0) the fluid fields {/, S} obey the complete 
set of hydrodynamic equations represented respectively 
by 



£1 
Dt 
dS_ 
'dt ' 



/V-V = 0, 



1 
2m 



|V5|- 



(17) 
(18) 



These are denoted as quantum hydrodynamic equations. 
The first one is manifestly the continuity equation for 
the quantum probability density J{v,t). Instead, the 
second one is the Hamilton- Jacobi equation for the 
quantum phase-function S. Moreover, Uqm is the so- 
called quantum potential |4| related to the interac- 
tion potential U by means of the equation Uqm = 



V2 In/ 



Vln/I 



U. Since by assumption 



U is single- valued in fi, it follows that equations H17() 
and (|18|l must also be single- valued. Nevertheless, both 
^ and Uqm are not unique since they are determined 
up to an arbitrary real smooth function z(t) since they 
are invariant with respect to the gauge transformation 



5(r, t) ^ S'ir, t) = 5(r, t) + i / dt'z{t'), 
U{r,t) ^ U'{v,t) = U{v,t) + z{t), 
f{v,t)~>f{r,t)^f{r,t). 



(19) 



S{rs,t) = Sn]{rs,t) mod{2TT) 



lim V(r,i) = V^(r5,i), 



(23) 



(24) 



where Sw{^s,t) and Vu,{v5,t) are suitably smooth real 
functions and fwi^s^t) is specified either by Eq.Q or 
(|10|l . depending on the definition of the domain $7. 

Equations H17() . H18|I . together with the initial condi- 
tions |j2ni),l|2lj and the boundary conditions 
define the quantum hydrodynamic initial-boundary prob- 
lem {QHE problem). 



III. GAUGE-INVARIANT FORM OF THE 
HYDRODYNAMIC EQUATIONS 

The gauge function z{t) can be eliminated by tak- 
ing the gradient of Ea. (|18|l term by term. The result- 
ing gauge-independent equations for the quantum fiuid 
fields {/, V} , again valid in the open domain SI x /, are 
provided by the gauge-invariant quantum hydrodynamic 
equations, which are defined by the continuity equation 
(fTTjl and by 

^V(r, t) = -¥ = --VUqm- (25) 
Ut m m 

As a consequence, by Eas. (|17|) and (|25|l can be viewed as 
the hydrodynamic equations of a compressible fluid. On 
the other hand, Ea. (|25|l implies 



d_ 

dt 



V X V(r, i) + V X (V • VV(r, t)) ^ 0, 



where V ■ VV = -Vx (V x V) - VV^ and 



V X (V- VV(r,<)) = 
= • VV -f V ■ 



-V X [V X 



(26) 



(27) 
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where ^ = V x V is the vorticity vector. Therefore, if 
we impose in the whole domain the initial condition 
^(r,to) = it results 

e(r,t) = (28) 

for all t G I. Notice that Ea. (|28|l is not in contradiction 
with the condition of quantization for the velocity circu- 
lation K [see Ea. (|16|l ] since the phase function S{r, t) re- 
sults generally non-single-valued. As a consequence, the 
vector field V(r, t) is necessarily of the form (|15() . Hence, 
the fluid described by the fluid fields {/, V =^VS'(r, t)} 
is necessarily vorticity-free, while at the same time its 
velocity circulation is non- vanishing [see Ea. (|16() ]. This 
equation is know as the so-called quantum Newton equa- 
tion '21 (or quantum Navier-Stokes equation 7\). The 
initial-boundary conditions for these equations are de- 
fined again by 1201), and which imply in 
particular for V(r, to) the initial condition: 

V{r,to) = -S/So{r,to). (29) 
m 

Equations (|17() . H25() ]. together with the initial con- 
ditions H2() |I . H21|I and the boundary conditions l|22l) - 
(I24|l . define the gauge-invariant quantum hydrodynamic 
initial-boundary problem (GI-QHE problem). 

In summary, by construction it follows that: 

a) the QHE problem is equivalent to the SE problem, 
namely {/, S} is a solution of the first problem if and 
only if ^(r, t) is a solution of the second one; as a conse- 
quence the solution {/, S mod(27r)} of the QHE problem 
is unique; 

b) if {/, S mod{2Tr)} is a solution of the QHE problem 
then {/, V =iV5'(r, i)} is necessarily a solution of the 
GI-QHE problem; 

c) vice versa, a solution {/, V} of the GI-QHE prob- 
lem, determines uniquely {/, S} up to an arbitrary gauge 
transformation of the form (|19|l . 

Finally, it is worthwhile to mention that, in principle, 
it is also possible to introduce sets of "reduced" hydro- 
dynamic equations, defined in appropriate subspaces of 
the A^— body configuration space, in particular the one- 
particle subspaces fli (for i = l,N). The latter can be 
obtained adopting for the A^— body quantum system the 
one-particle reduced representation described in the Ap- 
pendix A [see Ea. (|S^ ]. Manifestly, these reduced de- 
scriptions are not equivalent to the full A'— body descrip- 
tion. For example, the A^— body system can be consid- 
ered as formed by A^ 1-body subsystem, one for each 
particle {j — 1, N). For each 1-body subsystem it is pos- 
sible to introduce a set of reduced quantum fluid fields 
{fjj^j — m^J^ji^j'^)} ' both defined in the set ilj x / 
and uniquely associated to the one-particle wave function 
-i/jj (r^- , i) by means of Eqs. ^ and (|15|l . It is immedi- 
ate to prove that the fluid fields {/j,Vj} for j = 1,A^ 



obey a set of fluid equations formally analogous to Eqs. 
(|17|l . (|25() . to be denoted as reduced hydrodynamic equa- 
tions, which can be viewed as describing the dynamics of 

an immiscible fluid mixture. 

IV. HEISENBERG THEOREM AND THE 
CONCEPT OF QUANTUM TEMPERATURE 

The set fluid equations Eqs. ifTTjl and for the 
quantum fluid fields {/, V} provide a complete descrip- 
tion of quantum systems. This means, in particular, 
that no other independent observable or fluid field is re- 
quired to describe the quantum state. However, for the 
subsequent analysis it is useful to introduce the concepts 
of quantum directional temperatures and quantum tem- 
perature, which can be defined by analogy with classical 
statistical mechanics and interpreted as additional quan- 
tum fluid fields. We intend to show that the definitions of 
these observables emerge from Heisenberg theorem. We 
recall that this is realized by means of the (Heisenberg) 
inequalities (holding for z = 1, 2, 3) 

((Ar,)')((ApO')>f , (30) 

or 

An Aft > ^. (31) 

Here the notation is standard. Thus, Ar^ = 

/ 9\ — / 9\ 1/2 

^(Ar,)') ,Aft - {{^P^) ) (for i = 1,2,3) are 
the quantum standard deviations for position and linear 
momentum, ^(Ar^)^^ , ^(Aft)^^ are the corresponding 
average quadratic quantum fluctuations, while Ar^ , Aft 
are respectively the quantum position and momentum 
fluctuations Ar^ = r.;— (r^) , Aft = pi — (pi) . Finally 
(Q) = (-0 I (3-0) = J^drf{r,t)Q{r,t) denotes the expec- 
tation value of a generic dynamical variable Q. As usual 
we identify the quantum linear momentum p with the 
linear differential operator 

p = -ihV (32) 

which acts on the functional space . It follows (p) = 
(-0 I ptp) = /o drfV = (P) , where 

P = mV (33) 

is the fluid momentum, while the expectation value of 
(Aft)^, upon integration on the set fl, reads 

((Aft)^) -^Jdrf id, \nff + {{d,Sf) - {d,Sf . 
n 

(34) 
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As is well-known Heisenberg theorem follows by invok- 
ing the identity 



/(r,t) = l, (35) 




> 1. 



Hence, Schwartz inequality delivers: 



d 



(Ar^y) / drf 



1 2 



— ln/(r,v,t) 
on 



> 1, 



(36) 



(37) 



where by definition ({Ari)'^\ = J drf [r^ — 



A peculiar aspect of the Heisenberg inequality 1)30(1 is 
that it can also be written in terms of the relative fluc- 
tuations Ad) Pi = Pi - Pi, which are defined with re- 
spect to the components of fluid momentum Pi — mVi 
(for i = 1,2,3), instead of the corresponding expectation 
values {Pi). The analogy is based on the fact that the 
concept of kinetic temperature, and the related one of 
directional kinetic temperatures [see Ea. (|48|l in the next 
Section] , is defined in terms of fluctuations with respect 
to the local fluid velocity V(r,t), instead its expectation 
value (V(r,<)} . In fact, it is immediate to prove that, 
by definition of the quantum linear momentum (|32|l . the 
following identity holds (for i = 1, 2, 3) 



(Ap^ 



A(2) 



Pt 



(38) 



where /(Ad)pi)^\ and /(A(2)p,;)^\ read respectively 



A(l)r 



4 



drf id.lnfY 



(39) 



^(A(2)p,)'^ ^ (id,Sf)-{d,S)\ (40) 

and hence can be interpreted as the average quadratic 
momentum fluctuations carried respectively by the quan- 
tum probability density / and the phase function S. 
The average quadratic fluctuations ^(A^^^p^)^^ afford a 

straightforward interpretation in terms of classical statis- 
tical mechanics, by introducing the notions of quantum 
directional temperature TQM,i{t) and quantum tempera- 
ture TQM{t), defined respectively: 

mTQM.(i) = ^(a(S«)'^, (41) 



TQM{t) - 2 XI ^QM,i{t), 

1=1,2,3 



(42) 



which can be viewed as constitutive equations for 
TQM.zit) and TQM{t). We notice that if = R^^ and 
'0(r, t) is dynamically consistent j^^, necessarily it must 
result TQM{t) > TQM,i{t) > in 57. In the remainder we 
shall assume, however, that it results TQM{t) > for all 
€ I also in the case in which is a bounded set. As a 
consequence the following modified Heisenberg inequality 
holds 



(Ar, 



mrQM,»(i) + /(A(2)p,) 



>V (43) 



A particular case is the one in which the condition of 
(quantum temperature) isotropy 



Tqm,! [t) — Tqm {t) 



(44) 



holds identically for i = 1, 2, 3. In the following sections 
for greater generality we shall require, however, 



Tqm Ait) ^ Tqm, jit), 



(45) 



for i j (with i,j — 1,2,3). In fact, generally for ar- 
bitrary quantum systems Eg . ((44(1 cannot be assumed to 
hold. 



V. PROBLEM: THE SEARCH OF AN INVERSE 
KINETIC THEORY FOR NRQM 

A. Basic assumptions, the kinetic correspondence 
principle 

The form of the quantum fiuid equations 1(17(1 and 1(25(1 
suggests that they can be obtained as moment equations 
of a continuous inverse kinetic theory, analogous to that 
developed recently for the incompressible Navier-Stokes 
equation [see ET]. 

Let us now pose the problem of searching an inverse 
kinetic theory for the Schrodinger equation, i.e., a kinetic 
theory yielding exactly, by means of suitable moment 
equations, the quantum hydrodynamic equations. The 
theory must hold for arbitrary (and suitably smooth) ini- 
tial and boundary conditions both for the wavefunction 
and the kinetic probability density. In particular, this 
involves also the search of a possible underlying classical 
dynamical system which determines uniquely the time- 
evolution of the quantum system. 

In the sequel let us consider, without loss of gener- 
ality, the case of one-body quantum systems; the the- 
ory here developed is applicable, in fact, with minor 
changes also for systems with N > 1 particles. For def- 
initeness, let us assume that the quantum fluid fields 
{/(r, t), V, TgM.i, for z = 1, 2, 3, } are respectively solu- 
tions of the GI-QHE problem [i.e., Eqs. 
((21(1 . I(22() - ((24() ] and imposing also the constitutive equa- 
tion ((41() . To restrict the class of possible kinetic models. 
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following the approach of ET, let us introduce a probabil- 
ity density g{x,t), with x = (r, v) , defined in the phase 
space r = n X U (where U = M'^^) and assume that 
it belong to the functional class {g{x,t)} of real func- 
tions which satisfy the following properties 1-4 (denoted 
together as Assumptions #1); more precisely, it is as- 
sumed that g(x, t) 

1. is non-negative and continuous in F x /, in partic- 
ular, is strictly positive and of class C^*^'^-*, with 
fc > 3, in r X /; 

2. y{r,t) G X / admits the velocity moments 
Mx[g] = Jijdv^g. with X(r,v,i) = 1,v,m| (for 
i = I, 2, 3), uu, uu^, In g, where v, u = v — V are 
respectively the kinetic and the relative kinetic ve- 
locities and Ui = Vi — Vi (for i = 1,2,3) are the 
orthogonal Cartesian components of u defined with 
respect to an arbitrary inertial reference frame; 

3. admits the velocity moments X(r,v,i) — l,v,u| 
(for i ~ 1,2,3). The latter are prescribed by im- 
posing a suitable set of constraint equations, to be 
denoted as kinetic correspondence principle, which 
relate the quantum fluid fields and the correspond- 
ing kinetic moments. For this purpose the following 
equations are assumed to hold identically, respec- 
tively in X /, 



(46) 
(47) 



/ = Mi[.g] = / rfv.g(r,v,0, 
Ju 

M2 [g] = / dvvg(r,v,i) = 

= V(r,t), 

and in / for i = 1, 2, 3 

T,{t) = A/3, [g] = 
= / dvmu-5(r,v,i) = TQM,t{t) > 0. 



(48) 



Consistently with H45(l . we shall generally consider 
Ti ^ Tj for i ^ j (with i,j = 1,2,3). Here 
the moments Ti = M^i[g] (for i — 1,2,3) and 
T = {T1+T2 + T3) /3 are denoted respectively the 
( quantum ) kinetic directional temperatures and the 
(quantum) kinetic temperature. 

Finally, let us impose an appropriate regularity 
condition for the fluid fields and the quantum force 
F(r,t). In particular, besides imposing that the 
fluid fields {/, V,Tqa/,i, for i = 1, 2, 3, } are solu- 
tions of the GI-QHE problem and satisfy the con- 
stitutive equations H41|) . let us require that they 
belong to the functional settings: 



/,V,r,(r,t)GC('=^2)(^^ j)^ 
T,(r,i),/(r,<) > 0, 

/,v,r,(r,t) ec(")(f^ X/), 

F(r,t) G C('='2)(f7 X /), 



(49) 



with fc > 3. 

The constraint provided by Ea. (|48|l implies that the 
kinetic directional temperatures Ti (for i = 1, 2, 3) are as- 
sumed position-independent. This assumption, although 
consistent with the definition of the quantum directional 
temperatures given above [see Ea. (|41|l ]. may in principle 
be avoided (see related discussion in Appendix B and at 
the end of Sec. 6). 

Furthermore, let us require that in the open set F = 
Q.xU the probability density ^(r, v,f) satisfies a Vlasov- 
type kinetic equation of the form [Assumption ^2) 



Lg{^,t) = 



(50) 



[inverse kinetic equation), where L is the Vlasov stream- 
ing operator 



9 d _ d d 



(51) 



and^ = ■|r+v-V + — - ir- is the Lagrangian derivative. 

at ot m d\ ^ ^ 

Here x = (r, v) and X = |v,^K} , where K(x,i), to be 
denoted as mean field force, is a suitably smooth real 
vector field. 

Since the correspondence principle defined by 
Eqs.|gni>,(EZI) and gHl must hold identically in the 
open set x / it follows that the moment equations 
for Mi[g],M2[g] and A/a^ [g] , for i 1,2,3, must 
necessarily coincide identically, respectively, with the 
quantum hydrodynamic equations H17|) and (|25|l and the 
constitutive equation for the directional temperature 
(|41(l . In addition, it is obvious these moment equa- 
tions must hold also for arbitrary quantum fluid fields 
satisfying assumptions H49|) . Such implications will be 
discussed in detail below (see THM. 1 and 2). 



B. The phase-space Schrodinger dynamical system 

Let us remark that the kinetic equation (|5U|I deter- 
mines uniquely the time evolution of the kinetic distri- 
bution function g(x,<) in the whole extended phase-space 
Fx/ and consequently prescribes uniquely also the quan- 
tum fluid fields in the set x /. In fact, it can also be 
cast, respectively, in the equivalent Lagrangian and inte- 
gral forms 



^.(x(0,t) 



-5(x(i),i) 



d K(x(<),t) 



9v(t) 



TO 



(52) 
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J{^{t),t)g{x{t),t) ^ g{Xo,to) 



(53) 



where g{xo,to) is the initial kinetic distribution function 
and the curves {x(t) ,t G 1} define suitable phase-space 
Lagrangian trajectories. Moreover, the map 

7(x„,t„) : Xo x(<) = x(Xo, to, t) (54) 
is the flow generated by the initial- value problem 



-x= X(x,.) 

■^(^o) — Xq, 



and 



J{x{t),t) = exp i y df 



(55) 
(56) 

(57) 



is its Jacobian 

Here we shall prove that, due to continuity of the ki- 
netic distribution function and of the Schrodinger dy- 
namical system, this equation holds identically in the 
closure Fx/, except possibly in the nodes r„ G 6fl, 
where /(r„,i) — 0. However, by construction, the limit 

lim J{K{t),t)g{x{t),t)^g{xo,to) (58) 

exists for all G including all nodes. Instead, one 
can prove that the limit limrfj-j^^s Ji^it),t) does not 
exist if is a node. In particular, by suitable def- 
inition of the vector field K(x,i), we shall prove that 
x(Xo, to, t) results suitably regular so that the set of maps 
(I54|l generates a (generally non-conservative) classical dy- 
namical system {7(x .to)^^ to be denoted as phase-space 
Schrodinger dynamical system. 

These include in particular - besides Assumptions #1 
and #2 - the hypothesis that the kinetic distribution 
function g{:x.o,to) (and its initial condition g(:Ko,to)) is 
suitably smooth in the whole set Fx/. For exam- 
ple, we let us require that X(x,i) e C('^'^'(F x /) and 
gixo,to) e C^'^'^) (r X /), with A: > 3. Then it follows that 
x(Xo,io,t) is a diffeomorphism. of class C^'^'^'^)(Fx/x/) 
and g{x{t),t) e C^'^^^^T x /), with A: > 3 while, so that 
the moments {/, V,Ti, i = 1,2, 3} are necessarily of class 
C('='2)(f7 X /) {Assumption #3). 

Manifestly the Schrodinger dynamical system, if it ex- 
ists, provides a deterministic description of quantum me- 
chanics since it advances in time both the kinetic prob- 
ability density and the quantum fluid flelds {/(r, i), V}. 
Thus, a fundamental issue is the question of the existence 
of a vector field K(r, v,t) which satisfies the minimal as- 
sumptions indicated above. To be more specific, how- 
ever, it is convenient to further specify the mathematical 
model imposing an appropriate set of assumptions (de- 
noted together as Assumptions #4), which include the 
following ones [ET]: 



1. the kinetic equation admits local Maxwellian kinetic 
equilibria for arbitrary kinetic moments and quan- 
tum fluid fields which satisfy the correspondence 
principle i46}) , iT^ , i4S\ ) and the regularity require- 
ments i4y\) ; 

2. The vector field X(x,t) is prescribed in such a way 
that it depends, besides on x, only on the fluid fields 
and suitable differential operators acting on them; 

3. the kinetic distribution function satisfies appropri- 
ate boundary conditions for the kinetic distribution 
function; 

4. the Heisenberg theorem is satisfied identically. 

To complete the specification of the inverse kinetic 
equation, however, it must be supplemented with suitable 
initial and boundary conditions. In particular since the 
kinetic correspondence principle must hold both on the 
boundary 5Vt and at the initial time to, this means that 
on the boundary 5T (which includes 5^1) bounce-back 
boundary conditions are imposed on the kinetic distribu- 
tion function 1331 . 



C. Bounce-back boundary conditions 

We intend to define boundary conditions for the ki- 
netic distribution function g(x,i) which are consistent 
with the Dirichlet boundary condition defined on the 
boundary 5^1 for the quantum fluid fields. Denoting 
Y-u,{Y5{t),t) = -^rs{t) the velocity of the point of the 
boundary determined by the vector rs(t) e SQ and as- 
suming |v — Vw| 7^ 0, let us introduce the unit vector 



and the variable 



C = [v - Vw] • b. 



(59) 



(60) 



Here a ~ ±1 and its sign is defined so that when b is 
a vector applied at the position it is always oriented 
inward with respect to the domain il. For sake of defi- 
niteness, we shall assume that Sfl is a piece-wise regular 
surface and that the vector b belongs to the open tangent 
cone to SQ in r„ which is oriented inward with respect 
to fl. Then, at an arbitrary position G Sfl, the sign of 
the variable ^ determines incoming and outgoing veloc- 
ity subdomains, defined, respectively, as the subdomains 
of velocity space U for which ^ < and ^ > 0. There- 
fore g{Ts,w,t)\^^f^ and g{vs,'v,t)\^^Q denote the incoming 
and outgoing kinetic distribution functions at position r^. 
These notations permit us to define properly the bound- 
ary conditions for the kinetic distribution function on 

For any boundary (5ri, position r„ e and for all 
nonvanishing vectors v — taking at the directions 
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specified above, we impose, respectively, for ^ > and 
^ < the boundary conditions for (7(r, v,t) {Assumption 
#5) 

5'(r5,v,t)|^>o = .g(r5,2Vw - v,i)|^^o, (61) 

5(r5,v,t)|^^o = g(r5,2Vw - v,i)|^>o; (62) 

to be denoted as bounce-back b.c. for g(Ti,t). It is im- 
mediate to prove that they are consistent with Dirichlet 
boundary conditions defined by Ea. (|24|l . In fact, there 
results: 

fV{rs,t)^ / dvvg{rs,v,t) = 
Ju 

= / dvvg{rs,v,t) + / dvvg{rs,v,t) = 
Ju Ju 



1 r^" 

-/ dvv{g(r5,v,i)+.g(r5,2Vw-v,<)}+ (63) 
^ Ju 

1 /■«>° 

+ - J dvv{g{rs,^r,t)+g{rs,2V^-v,t)}. 



lu 

Thanks to the identities 
i J dvvg{rs,2V^ -v,t) 



Vw-2V(r„i) 



/ dvvg(r5,v,i) = V(r„i)/ 
Ju 



f (64) 



(65) 



it follows that: 



yw{rn,t)^^ [ dvv.g(r„,v,i) = V(r„,t). (66) 
/ Ju 

We impose, furthermore, the Dirichlet boundary condi- 
tion [Assumption #6) 



dvg{rs,v,t) = f^{rs,t), 



(67) 



where (r^ , t) is the prescribed probability density on 
the boundary defined by (j221l- 



VI. CONSTRUCTION OF THE MEAN-FIELD 
FORCE K 

A. Case of the generalized Maxwellian distribution 

It is well known that the principle of entropy maxi- 
mization (PEM) 13 113, m 

can be invoked to deter- 
mine uniquely, at a prescribed initial time to G I, the ini- 
tial condition for the kinetic distribution function g{:x., t). 



This requires that the functional class {(7(x, t)} must 
be suitably specified; in particular PEM can be used 
to determine the initial condition g{'Si., to) which corre- 
sponds to the requirement that quantum fluid fields ( 
{/(r, i), V,rQM,i, for i = 1, 2, 3, } are uniquely prescribed 
at the initial time t = to- In such a case the kinetic 
correspondence principle determines uniquely the kinetic 
distribution. This reads necessarily (for t — to) 



5M(r,v,t) = f{r,t) 



1 



exp{-a;jXi} (68) 



(generalized Maxwellian distribution), where in the ex- 
ponential the sum is understood on repeated indexes 
{i — 1,2,3). Here the notation is analogous to ET, thus 
for 1 = 1,2, 3, vthi = \/2Ti(r,t)/m and Xj = Ui/vthi- 

Since the initial time to is arbitrary, it is natural to 
assume that gA/(i",v,t) results identically (V(x,t) e F x 
/) a particular solution of the inverse kinetic equation 
(|5n|l . In such a case, invoking Assumptions #l-#6 (and 
in particular the hypothesis that the kinetic directional 
temperatures can only depend on time), the mean-field 
force K(gM) results necessarily (Assumption #7) 



with 



K{gM) = Ko(5m) + Ki(5m), 



Ko(ffAf) = F(r,t)-f iy^;Le.e,- V/ 

1, 



(69) 



(70) 



F(r,t) + -Vp, 



Ki(g 



M 



mu • VV- 



— w,e, — InT,-. 
2 ' dt 



(71) 



where the sum is understood on repeated indexes and p — 
JT denotes the kinetic scalar pressure. Here, Ko(5m) 
and 'Ki{gM) have been distinguished for being, respec- 
tively, constant and velocity dependent. In particular, 
Ko(ffM) contains, besides the quantum force F(r,i), a 
corrective term (" pressure term" ) which depends explic- 
itly only the logarithmic gradient of /; instead Ki((;m) 
contains a " convective term" , proportional to VV and a 
contribution proportional to the logarithmic time deriva- 
tives of the directional temperatures. For the sake of 
reference, the more general case in which the kinetic di- 
rectional temperature are taken as spatially non-uniform 
is reported in Appendix B. 

Let us now examine the main implications which 
stem, in the particular case g = ^m, from positions 
(|5^ - (f7T|) and assumptions #l-#7. We first notice 
that the Schrodinger dynamical system generated by 
K(gM), defined as the solution of the initial value- 
problem H55|l . (|55| l. exists and is unique in the whole ex- 
tended phase-space Fx/. This property is manifestly 
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assured by assumption, thanks to to previous definition 
of K((7m) and the regularity properties of the quantum 
fluid fields (Assumption #1). In the same set the ki- 
netic distribution function (7A/(r, v, t) is, by construction, 
a particular solution, i.e., exists, is unique and results 
strictly positive. In fact, this property holds if and only 
if the quantum fluid fields are solutions of the GI-QHE 
problem and - in validity of Eas. (|69() - (|71() - if the kinetic 
directional temperatures are assumed to be only func- 
tions of time, i.e., Ti(t) (« = 1, 2, 3). In addition, by con- 
tinuity, the kinetic distribution function (7M(r, v,i) (and 
its moments) are uniquely defined also on the boundary 
set Sn and hence, in particular, in the nodes. Finally, it 
is immediate to prove that Ea. H5()(l is also an inverse ki- 
netic equation for the quantum hydrodynamic equations. 
Indeed its moment equations, evaluated with respect to 
the weight functions G(x,i) — l,v, coincide identically 
(in O X J) with the Eas. H17|l . (|25() . As a consequence, the 
following theorem holds for Maxwellian kinetic distribu- 
tions of the type (|68|l : 

THEOREM 1 Generalized Maxwellian solu- 
tion of the inverse kinetic equation 

Besides the validity of Assumptions #l-#7, let us re- 
quire that: 

1 ) the kinetic distribution function fulfills Eq. at 
least at the initial time to G /; 

2) the mean-field force K can depend functionally but 
not explicitly on qm] moreover KolgM) and Ki{gM) can 
only depend on {f{r,t),'V,Ti,for i = 1,2,3,}. 

Then it follows that: 

A) y (x,i) e r X / the generalized Maxwellian distribu- 
tion HbfAl is a particular solution of the inverse kinetic 
equation l^5(J\) if and only if the mean field force K has 
the form defined by Eas. ^69\) . j70| ) and ^71}) ; 

B) the Schrodinger dynamical system exists, is unique 
and is continuous inTxI, except in the nodes, i.e., when 
r{t) = r„, with f{r{t),t) = 0. Moreover, it is C^'^^'^^TxI) 
with fc > 3. This result holds for arbitrary quantum dy- 
namical systems and both for isotropic and non-isotropic 
quantum directional temperatures; 

C) the Jacobian of the phase- flow Xq ^ x = x(t) — 
x(Xo,io,t), generated by the initial value problem 
\55]) .\5b ]) . is defined V(x,i) £ F x /, except in the nodes. 
There results in such cases 

Ki{Yo,to)f{'c{t),t)eii-p{-x^{t)x^{t)] 7^ (72) 

so that J{:x.{t),t) reads : 

j( (,^ ,N ^ Ki{t)f{{ro,to))exp{-x,oX,o} 

Kl{to)f{r{t),t)exp{-x^{t)x^{t)}' ^ ' 

where Ki = iTiT2T3)^^^ , Xoi = Uoi/ vth,i{to), Uo = 
Vo-V(r„,io), x.,{t) = Ui{t)/vth.i{to) and u(i) = 
v(t)-V(r(t),i); 

D) the limit 

lim J(x(i),t)5M(x(t),t) (74) 

r(t)-r„ 



exists and is unique; 

E) V(x, t) G F X /, the velocity-moment equations of 
the inverse kinetic equation \5U\) evaluated for the weight 
functions G{x,t) = 1, v and for g = gM coincide with the 
fluid equations Eqs. il7}) , ^25\) : 

F) V(x, f) G F X /, the moment equations for the direc- 
tional kinetic temperatures are satisfied identically. 

PROOF: 

If the mean- field force K is assumed of the form H69|l . 
(|70|1 and (|71|l the proof of A follows by straightforward 
algebra from Ea. H50|l . or in an equivalent way from the 
integral equation Ea. (|53|l . Fhis furnishes also the proof 
of Eq. (|73|l , which follows invoking Liouville theorem H57|) 
(C). Moreover, since all moments /, V, Ti {i = 1,2,3) 
and the quantum force F(r,i) by assumption belong to 
the functional setting H49I) , the proof of B is an immediate 
consequence of the fundamental (existence and unique- 
ness) theorem for ordinary differential equations. In par- 
ticular, it is obvious that the solution of the Schrodinger 
dynamical system is not defined in the nodes since the 
mean- field force K is not defined for f — 0. Regarding 
C, the only possible singular behavior of J(x(t),i) can 
occur either in the nodes r„ G (SO, i.e., if at some time 
ti G / there results /(r„,ti) = 0, or if least one of the 
directional temperatures Ti vanishes. The second possi- 
bility is excluded by assumption. Hence J(x(t),i) is not 
defined at (x(i),t) only if r{t) = r„ is a node. Neverthe- 
less, it is obvious that the limit (|74|l exists (D) and more- 
over that 5a/ (x, t) is by construction non negative and is 
manifestly defined everywhere in Q x I. Finally the mo- 
ment equations satisfied by the inverse kinetic equation 
for g = gAf(r,v,i) are straightforward. It is immediate 
to prove that the first two moments G = 1 , v coincide 
respectively with the quantum hydrodynamic equations 
((131,(123, while G = uf (i = 1, 2, 3) yield moment equa- 
tions for the directional temperature which are satisfied 
identically (E,F). 



B. Case of the non-Maxwellian distributions 

As in Ref. ET the inverse kinetic theory can be formu- 
lated for non-Maxwellian kinetic distribution functions 
too. In fact, although PEM determines uniquely H()8|l . 
different initial conditions are conceivable in which the 
kinetic distribution function (7(r,v,io) = 9o(j,'v) is an 
arbitrary function, different from (|68(l . but results oth- 
erwise, by assumption, suitably smooth, strictly positive 
and summable in F. In analogy with ET, a unique def- 
inition of K((7) which yields the correct fluid equations 
and satisfies also the constitutive equations can readily 
be obtained also in this case. There results (Assumption 
#7b) 



K(g) = Ko(g) + Ki(5), 



(75) 
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Ko(5) = F(r,t) 



Ki(.g) = TOU • VV- 



m ^ 

- — UiBi 

2 



^1 T 



(76) 



(77) 



where again the sum is understood on repeated indexes. 
Equations (|75() , H77I) hold if the kinetic directional temper- 
atures are assumed to be only functions of time [Ti{t), 
for i = 1,2,3]. In this case the corrective term in the 
mean- field force Ko((7) depends on the tensor pressure 
n, instead of the scalar pressure p, while Ki((/) contains 
an additional term depending on the relative heat fluxes 
Qi {i — 1,2,3). The kinetic moments n and Qi (for 
i = 1,2, 3), assumed to exist, are respectively 



(78) 
(79) 



n = y dvuu^. 



The mean- field force defined by Eqs.((71I),l(7ni,(IZ3l is 
manifestly consistent with the previous definition when 
9 = ffM, since in such a case there results n = fTl ~pl 
and Qi= (i = 1,2,3). ~ 

It is immediate to prove that also in this case the 
Schrodinger dynamical system generated by ^.{g) exists 
and is unique in the whole extended phase-space Fx/. 
In the same set the kinetic distribution function g{r, v, t) 
exists, is unique, results strictly positive and by conti- 
nuity is uniquely defined also on the boundary set 5il. 
Finally, by construction the moments of the inverse ki- 
netic equation (|5n|l , evaluated with respect to the weight 
functions G(x,i) — l,v, coincide identically (in fl x I) 
with the quantum hydrodynamic equation Eas. (|17(l . l25|l . 
As a consequence, in the case of non-Maxwellian kinetic 
distributions the following theorem holds: 

THEOREM 2 Non-Maxwellian solutions of 
the inverse kinetic equation 

Besides the validity of Assumptions #l-#7, with #7b 
replacing #7, let us require that: 

1 ) the kinetic distribution Junction is a smooth func- 
tion of class C'^^''^^{T X /), with k >2, which is suitably 
summable in F; 

2) the mean-field force K can depend functionally 
but not explicitly on g; moreover Ko{g) and Ki{g) de- 
pend, besides {f{r,t),'V,Ti,fori — 1,2,3,}, at most on 
the higher-order moments 11 and (i = 1,2,3). 

Then it follows that: 

A) if at a prescribed time t G I and Vx e F the kinetic 
distribution function g{x-,t) is of the form then it 
results g(^,t) — gM{x.,t), V(x,i) G F x /; 

B) V(x, <) G F X / and for arbitrary smooth g{K-,t) the 
velocity-moment equations of the inverse kinetic equation 



\5(]\) evaluated for the weight functions G{'K,t) = l,v co- 
incide with the fluid equations Eqs. il7}) , ^25\) : 

C) V(x, t) G F X /, the moment equations for the direc- 
tional kinetic temperatures are satisfied identically; 

D) the Jacobian of the phase-flow Xo — > x = x(t) — 
x{^o,to,t), generated by the initial value problem 
<55|) .i)56 ]) . is defined V(x,<) G F x / reads in this case: 



J(x(t),i) = 



K,{r{t),t)f{ro,to) 
K^{vo,to)f{v{t),t) 



(80) 



expjy dt'G{x{t'),t') 

where Ki = {TiT2T3)^^^ and 

G(x(t),t) = iu.V/+i^V.Q- (81) 

E) if g{'K,to) is strictly positive V(x) G F, then for all 
V(x, i) G F X /, g{x,t) is also strictly positive; 

F) the Schrodinger dynamical system exists and is 
unique and is continuous in T x I, except in the nodes, 
i.e., when r(t) — r„, with f{r{t),t) = 0. Moreover it is 
(7(fe,2)j^p ^ j-j with k > 3. This result holds for arbitrary 
quantum dynamical systems and both for isotropic and 
non-isotropic quantum directional temperatures. 

PROOF: 

The proof of A is immediate thanks to THM.l, since 
by construction Ti{g) ~ K.{gM). The proof of B and C 
follows by direct inspection. In particular, again the di- 
rectional temperatures Ti (i — 1,2,3) are completely ar- 
bitrary, but suitably smooth, time-dependent functions; 
hence they can be uniquely determined according to 
Eq. (|^ . Eq. l|5n|) in D follows after straightforward alge- 
bra invoking Liouville theorem (|57|l . In particular, once 
again it follows that J(x(t),t) is not defined if r{t) [or r„] 
coincide with a node r^, G Sfl where /(r„,i) = 0. As a 
consequence V(x, t) G F x /, J(x(t),t) > 0. This proves 
also propositions E and F. Finally, to reach the proof of 
uniqueness of the mean-field force T^{g) it is sufficient to 
prescribe, besides the assumption that its does not de- 
pend explicitly on g, also an appropriate dependence in 
terms of higher-order kinetic moments. 

An important consequence of the previous theorems is 
the uniqueness of the mean- field force K((7) and hence of 
the inverse kinetic equation. 

THEOREM 3 Uniqueness theorem for K{g) 

In validity of THM.l and 2 the mean field force, for 
which A-E hold, K((7) is unique. For g = gM it has the 
form defined by Eqs. 1^6 9\) . j70| ) and j7J| ); for g ^ gM the 
mean field-force has the form defined by Eqs. j75| ), j76p 
and {73). 

PROOF 

For example, let us prove that in validity of assumption 
3) of THM.l the mean-field force K((7m) is unique. Let 
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us assume that there is a nonvanishing vector field AK 
such that K' = K+AK is also an admissible mean-field 
force. Hence it must be 



d_ 
dv 



(AK^m) = 0. 



(82) 



This means that AKg^ has necessarily the form 



dA dB 



(83) 



where A and B are suitably smooth real scalar fields. If 
both A and B are independent of gM then it must be 
AK = 0. Instead, for example, let us assume that only 
A depends on qm- Letting 



1 dA{gM,r,v,t) 



9M 



9v 



(84) 



if follows that A must depend on {gM too. Hence, in 
order that AK results independent of gM it must vanish 
identically. This proves that it must be AK = and 
hence K is unique. 

Let us briefly comment on these results. 

First, we stress that Schrodinger dynamical system 
[defined by the initial- value problem H55|) . l|56|l ] uniquely 
generates the time evolution of the kinetic distribution 
function g and hence of its moments. The result holds 
not only for generalized Maxwellian distributions H68|l 
but also for arbitrary, but otherwise suitably smooth, 
kinetic probability densities (THM.2). An interesting as- 
pect concerns, in particular, the behavior of the kinetic 
distribution function and of the Schrodinger dynamical 
system in the nodes, i.e., the points of the boundary Sil 
in which the probability density /(r,t) vanishes. Their 
occurrence is usually associated to the possible singu- 
lar behavior of the wavefunction, i.e., its non-uniqueness 
(for a review see Therefore, the question arises 

whether the existence of these roots is reflected in any 
way in the inverse kinetic approach (and therefore on 
the quantum system) . It is manifest that in these points 
the dynamical system is not defined and hence the Ja- 
cobian of the Schrodinger dynamical system [given by 
(17!^ or more generally by Ea. lp^ ] vanishes or is not de- 
fined. However, in these points the kinetic probability 
density can still be uniquely defined, thanks to conti- 
nuity. Consequently, the kinetic distribution function g 
and its moments, and in particular those corresponding 
to the relevant quantum fluid fields, are defined without 
exceptions in the whole domain f2 x /. As a consequence, 
it follows that no singularity appears, i.e., the quantum 
fluid fields (and hence the quantum wavefunction) are 
unique in fl x I. 



C. The kinetic representation of Heisenberg 
inequalities 

It is interesting to examine the role of the kinetic di- 
rectional temperatures (i = 1,2,3) in the framework 
of the inverse kinetic approach and the consequent inter- 
pretation of the Heisenberg theorem. For this purpose, 
let us analyze the implications due to the kinetic corre- 
spondence principle. It is immediate to prove that the 
Heisenberg inequalities (|31l) can be represented, in terms 
of the kinetic standard deviations for position and linear 
momentum Ar*^™ , Ap^^" , in the form: 



Arf"Apf" > 



(85) 



^km g^j^^ Ap^™ are defined in terms of 
appropriate phase-space averages. For this pur- 
pose, introducing the phase-space average {{A)) — 
J dxg{r,v,t)A{r,v,t),w'heTe A{r,v,t) is an arbitrary 
r 

summable phase-space function, let us pose respectively 

2\ \ 1/2 



r— - u(Ar,)'))"", Apf" = ( ( (Apf 



(for i = 1,2,3), where (({Ar,) )) and Ap: 



. —kill 

Ap, 



are denoted as average quadratic kinetic 



fluctuations. Here Ar = r— (r) is the position fluctua- 
tion, while Ap*^™ = p''™— to(V) and p*^™ = mv are 
respectively the kinetic momentum fluctuation, and the 
kinetic momentum. It is immediate to prove that (|85|l are 
equivalent to (|31|l . In fact, in analogy with Ea. (|38l) . one 
finds that the average quadratic momentum fluctuation 



A —kin 

APi 



can be written in the form 



kin 



(A(2)p,) \ (86) 



[i = 1,2,3), where Ap**" denotes Ap 
pfe»n_^Y(r,i) and ^(A^^)^^)^^ is given by Eq.gUJl. 
Moreover, by definition, it follows that the momentum 



fiuctuation 



reads 



kin \ 



mTi{t), 



(87) 



where and Ti{t) (for i — 1,2,3) are the directional ki- 
netic temperatures. As a consequence, the constraint 
Ti{t) = TQM,i{t), set (for i = 1,2,3) by Eq.lgHl of the 
correspondence principle, implies 



Ad) 



Pi 



= (((Ap? 



(88) 



It follows that the Heisenberg inequalities can be inter- 
preted in terms of kinetic fiuctuations, i.e., as constraints 
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between (J{^(Ar,)^)) and Apf" = ^{^(Apf") , 

whereby the quantum observable p — —ihX/ and its av- 
erage quadratic quantum fluctuation are replaced by by 
the kinetic momentum p'^*" — mv and the correspond- 
ing average quadratic kinetic fluctuation. This result is a 
direct consequence of the assumption set by Eq. for 
the directional kinetic temperatures Tj (i = 1, 2, 3). 

D. Generalizations: non- uniqueness 

It is obvious that the present results can be general- 
ized in several ways. In particular, the definition of the 
kinetic directional temperatures remains in principle ar- 
bitrary since they do not enter explicitly the quantum 
hydrodynamic equations. Thus, for example, it is possi- 
ble to require that the functions Ti {i — 1, 2, 3) are also 
position-dependent (see Appendix B). 

Due to the arbitrariness of the kinetic temperatures, 
it follows that there exist infinite equivalent realizations 
of the Schrodinger dynamical system and of the associ- 
ated Lagrangian trajectories {x(i), t e /} which yield the 
same quantum hydrodynamic equations. Therefore, the 
unique inverse kinetic theory here presented, which corre- 
sponds to a well-defined set of prescriptions and in par- 
ticular the assumption of spatially-constant directional 
temperatures, is simply one of the infinite possible math- 
ematical realizations. 

A side aspect concerns the so-called uniqueness prob- 
lem of the deterministic viewpoint of SQM 0, lilt IT^ , 
i.e., the Bohmian program, of reproducing the predic- 
tions of SQM within the framework of suitable determin- 
istic Lagrangian trajectories. In fact, the Schrodinger dy- 
namical system determined by the initial-value problem 
(|55|l . (156(1 yields in terms of the associated Lagrangian tra- 
jectories {x(i), t G 1} a deterministic description of SQM. 
Hence, it can also be viewed a phase-space generalization 
of Bohmian mechanics ^ lol - An implication of the 
present theory is that such a program has by no means 
a unique solution. In fact, there are infinite equivalent 
Lagrangian trajectories determined as solutions of the 
initial- value problem ((55|l . ((56(l . which differ only by the 
choice of the spatial-dependency assumed for the direc- 
tional temperatures (see Appendix B). 

VII. CONCLUSIONS 

Motivated by the analogy between hydrodynamic de- 
scription of SQM and classical fluid dynamics an inverse 
kinetic theory has been developed for the quantum hy- 
drodynamic equations. We have shown that, although 
in principle infinite solutions to this problem exist (in 
particular due to the indeterminacy in the kinetic di- 
rectional temperatures), the inverse kinetic theory can 



be uniquely determined, provided appropriate hypothe- 
ses are introduced. The results presented are relevant for 
the fluid description of quantum mechanics and a deeper 
understanding of the underlying statistical (in particular, 
kinetic) descriptions. 

For this purpose the full set of gauge-invariant quan- 
tum hydrodynamic equations, including the Heisenberg 
inequalities, have been related to the appropriate quan- 
tum fluid fields. As a result the notions of quantum 
temperature and quantum directional temperatures have 
been introduced. The present approach has the following 
main features: 

1 . the inverse kinetic equation 1)50(1 has been assumed 
to be a Vlasov-type kinetic equation; 

2. its solution, i.e., the kinetic distribution function 
has been required, in particular, to admit kinetic 
directional temperatures Ti (i = 1, 2, 3) which, con- 
sistent with the correspondence principle [defined 
by Eqs.(jin|),l|lZIl and ^^], depend only on time; 

3. the inverse kinetic theory holds for suitably smooth, 
but otherwise arbitrary quantum fiuid fields, which 
satisfy the quantum hydrodynamics equations with 
appropriate initial-boundary value conditions; 

4. by imposing a kinetic correspondence principle, i.e., 
by identifying the relevant quantum fiuid fields with 
appropriate kinetic moments, the quantum hydro- 
dynamic equations are satisfied identically, when 
are expressed in terms of the relevant moment equa- 
tions. 

5. theory is non-asymptotic, i.e., the quantum hydro- 
dynamic equations are satisfied exactly; 

6. the moment equations form a complete system of 
equations (closure condition). 

7. under suitable assumption, the inverse kinetic the- 
ory and the mean-field force which defines the 
streaming operator are unique. 

An interesting result of the theory, relevant for the 
mathematical investigation of the Schrodinger equation, 
concerns the discovery of the underlying dynamical sys- 
tem, i.e., the phase-space Schrodinger dynamical system. 
We have found that this can be identified with the non- 
conservative dynamical system advancing in time the ki- 
netic distribution function and generated by the kinetic 
equation itself. The evolution of the fluid fields is proven 
to be determined uniquely by this dynamical system. 
Formally the Schrodinger dynamical system can be in- 
terpreted as describing the dynamics of system of clas- 
sical "virtual" subquantum particles which interact with 
each other only by means the mean-field force K and are 
characterized by a dynamics which fulfills a suitable set 
of regularity assumptions. 
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To conclude, a further interesting feature of the present 
treatment is its adoption of the mean field force kinetic 
model. This permits, in principle, numerical implemen- 
tations by means of appropriate algorithms based on the 
new inverse kinetic theory [see related discussion in ET] . 

A side aspect concerns also the uniqueness and regular- 
ity of the quantum fluid fields and of the related quantum 
wavefunction. We have shown, in fact, that the kinetic 
distribution function is smooth in the whole extended 
phase space T x I, while is uniquely defined also in the 
nodes, i.e., the "singular" points of configuration space Q 
where the quantum probability density /(r,f) vanishes. 
As a consequence, the quantum fluid fields are necessarily 
unique in f2 x / and suitably smooth in O x /. 

Finally, have pointed out that the Heisenberg inequal- 
ities afford a simple statistical interpretation, which per- 
mits the representation of the quantum statistical fluc- 
tuations of the components of the linear momentum, 
and corresponding quantum directional temperatures in 
terms of statistical fluctuations of the kinetic momentum 
and of the kinetic directional temperatures. 
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prescribes the state of the j-th one-particle subsystem 
is defined by means of the integral 



C(,)(i)V',(rj,i) = L,ib^^Hr,t) EE ^j(r„i). 



where Lj is the integral operator 



dri..drN 



(90) 



(91) 



n 



Here Cj{t) (for j — l,N) are a real functions defined so 
that there results identically 



1 



cm 



dvi 



V^,(r„i) =1 (92) 



and fj — (r^ , t) \ are the associated one-particle prob- 
ability densities. The Schrodinger equation for ^j(rj,t) 
follows immediately from the A'^— body Schrodinger equa- 
tion (PJ. Let us assume for definiteness that the A^— body 
Hamiltonian takes the form 



k=l,N k,m=l 
Uokirk,t), 



k+ ? . Uk m I'm; t) + (93) 

k.m—l.N^kKm 



(94) 



where f7fc,„(rfc, r,„, i) and Uok{i^k,t) are respectively bi- 
nary and unary interaction potentials. Introducing the 
position 



APPENDIX A: REDUCED ONE-PARTICLE 
DESCRIPTION OF QUANTUM SYSTEMS 

In SQM the state of a system of N interacting particles 
is, by assumption, represented by its A'^— body wavefunc- 
tion ^p{r,t), with r = (ri,...,rAr) G Q, rj £ Qj (for 
j = 1, A^) and t & I. This is defined in the set fl x I, 
where O is the configuration space fl = Y[ > ^ 

j=l,N 

(for j = 1, A^) and / is an open subset of R. However, 
since the number of particles forming a quantum dynam- 
ical system is "a priori" arbitrary, also "reduced" quan- 
tum descriptions of an A^— body system are permitted. 
These descriptions, however, are not equivalent to the full 
A^— body description based on the A^— body wavefunc- 
tion. Thus, it is in principle possible to obtain a reduced 
description based, for example, on one-particle wavefunc- 
tions, whereby the A^— body system is represented by the 
reduced vector state 



ipBir,t) = (V'i(ri,t), ....ipNirN,t)), 



(89) 



instead of the single scalar A^— particle wavefunction 
ij}{r,t). The one-particle wavefunction 'ijjj{rj,t) - which 



= 1,N 

cu){t)Ui{rj,t)il;j{rj,t), 



J ^ C/fc„(rfc,r„„t)V'<'^^(r,t) = (95) 

k ,m— \ \ h<.'rn 



(96) 



it follows 

L.Hi, = C(,)(t) 4- [/i(r„t) -f L/oj} i^,{v,,t). (97) 

Hence, ipj{rj,t) obeys necessarily the one-particle 
Schrodinger equation 



d 



(98) 



where the index j (for j — l,N) identifies the particle 
subsystem (or species index) and 



- Hoj 



d 

Uoj - ifi-g^^'ac(^j) (99) 



is the j-th particle Hamiltonian. The reduced one- 
particle description of a A^— body quantum system is, 
therefore, obtained by means of the vector wavefunction 
■0(r, i) = {ipi, ....iPn) , ^j{Yj,t), for j = 1,N, being the 
one-particle wavefunctions which obey Ea. H98|) . 
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APPENDIX B: CASE OF 
POSITION-DEPENDENT DIRECTIONAL 
TEMPERATURES 

We notice that the correspondence principle H48|) can 
be modified by assuming instead Ti = 7i(r,t) {i — 1, 2, 3) 
and imposing, in place of Eq. H48|l , the constraint equation 

(T,(r,<)) =Tqm.W, (100) 
with general solution of the form 

T,(r,<) = fc(,)(r,t)(r,(r,t)}, (101) 
(fc(,)(r,t)) = l. (102) 

The functions fc(i)(r,i) (j = 1, 2, 3) which satisfy Eg. ((117^ 
are manifestly non-unique. In this case it is immediate to 
prove that for the generalized Maxwellian solution H()8|) 
the mean-field K{gM) is obtained by imposing (|^ . (|7n|) 
with 



Ki(5A/) = rriu- W+—Uiei—lnTi- 
-fe,e,<,. VlnT,(a:2-i)- (103) 

replacing Eq. H71|l . Here, again and the sum is under- 
stood on repeated indexes. The general case in which 
9 7^ 9m can be obtained immediately from Eas. (|103|l and 
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